Dynamics and steady states in excitable mobile agent systems 
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We study the spreading of excitations in 2D systems of mobile agents where the excitation is 
transmitted when a quiescent agent keeps contact with an excited one during a non-vanishing 
time. We show that the steady states strongly depend on the spatial agent dynamics. Moreover, 
the coupling between exposition time (u) and agent-agent contact rate (CR) becomes crucial to 
understand the excitation dynamics, which exhibits three regimes with CR: no excitation for low 
CR, an excited regime in which the number of quiescent agents (S) is inversely proportional to CR, 
and for high CR, a novel third regime, model dependent, where S scales with an exponent ^ — 1, 
with ^ being the scaling exponent of lj with CR. 

PACS numbers: 87.23.Cc, 87.19.Xx, 02.50.-r, 87.18.Bb 



The understanding of information propagation 
through a system of moving agents is crucial for many 
applications, ranging from chemical reactions to epi- 
demic spreading [ij. A particular example of such a 
process is the propagation of an excitation through a 
system of mobile agents. In an excitable system, agents 
typically have three states, quiescent, excited, and 
refractory, and pass through them in such a way that 
once excited become first refractory, and then quiescent 
again. A classical example of an excitable system is 
the forest-fire In this case, the immobility of the 
agents is well justified. Another example of excitable 
system is a disease in which agents undergo a cycle by 
which an infected agent becomes first recovered and then 
susceptible (SIRS dynamics). Though in the spreading 
of diseases the mobility of the agent has been for a long 
time ignored JsMjL, its importance has been recently 
recognized (aSM, ■ The first attempt to incorporate 
interactions among individuals was through the study of 
disease spreading on complex networks which typi cally 
represent social structures 13, 

The next step was to incorporate the actual motion 
of individuals across subpopulation structures where 
nodes represent containers of individuals such as urban 
areas or just communities 0, @, 0, H, 0. Despite 
the fact that complex networks can describe realistic 
social interactions or traffic networks, these models 
constitute frozen pictures where edges do not exhibit 
any dynamics. Gonzalez et al. in [l^ showed that 
the degree distribution and clustering coefficient of 
empirical social networks can be described by a systems 
of off-lattice mobile agents, provided the agent density is 
appropriately chosen. Interestingly, in this model agents 
are constantly moving and the underlying agent-agent 
contact network is highly dynamical. 




While agent motion may be central to informa- 
tion/excitation spreading, however, it poses questions 
that are not easy to answer. How does the spatial agent 



FIG. 1: (a) simulation snapshot. Number of agents A'^ = 1024. 
Colors indicate agent states: quiescent/susceptible (green), 
excited/infected (red) and refractory/recovered (blue). A 
short segment indicates the active direction of motion, (b) 
amplification of the area indicated in (a), (c) temporal evo- 
lution of the populations. 



dynamics affect the excitation spreading? In a SIRS- 
dynamics, is the agent-agent contact rate (CR) the key 
parameter that controls the the epidemic dynamics? 

In classical epidemiological mod els Id, \M as well as in 
previous mobile agent models [l^, llll. Il2l. 13 1, it is as- 
sumed that the mean number of neighbors (k) (directly 
associated to CR) and the effective infection probability 
A are independent. Consequently, for either a SIRS or 
SIS dynamics there are only two regimes with CR, resp. 
(k). For low CR the system remains in the epidemic ex- 
tinction regime, while for CR above a critical value the 
disease propagates such that the number of susceptible 
agents is inversely proportional to CR. 

In this paper we study for the first time an excitable 
mobile agent system where the excitation is transmitted 
by keeping physical contact for a finite time (see Fig. 
ID). We derive expressions for the threshold and steady 
states of the excitation dynamics, which highly depend on 



2 



the spatial agent dynamics, and show that the coupHng 
between exposition time (which directly affects A) and 
CR (resp. (k)) is crucial to understand the evolution of 
the excitation. In particular, we show that due to this 
coupling the excitation dynamics exhibits three regimes 
with the CR. These results applied to epidemic spreading 
suggests that, contrary to what classical epidemiological 
models state, the epidemic dynamics exhibits more than 
two regimes with CR. 

Spatial agent dynamics. — We assume agents are self- 
propelled disks which, in absence of interactions, move at 
constant speed in a box with periodic boundary condi- 
tion and change their direction of motion at Poissonian 
distributed times. The equation of motion of the i-th 
agent can be expressed in the following way: 

x.(t) = ^ + i^VC/(x,(t),x,(t)) (1) 



c c 



where F^/^ is the active velocity of the agent and ( de- 
notes a friction coefficient. As mentioned above, the 
active direction of motion follows a Poisson process by 
which an agent changes the direction of with a rate a 
while keeping constant the active speed v — \Fi \ /(^ (loj . 

As example of agent-agent interaction we have chosen 
a repulsive soft-core two-body potential that penalizes 
agent overlapping: 



t/(x,x') = 
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|x-x'|-^ - (2r) 



-0 



if |x-x'| < 2r 



if |x-x'| > 2r 



(2) 



where r is the radius of the agents, /3 is a constant, and 
7 is a function of v such that the maximum overlap- 
ping area between two agents is independent of v [l9|. 
The collision between two or more agents is a relatively 
slow process in which agents keep physical contact for a 
non-vanishing time. Notice that while the implemented 
dynamics is suitable to describe migration of some micro- 
organisms I20], it becomes unrealistic as a model for hu- 
man mobility. 

Excitation dynamics. — Without lack of generality we 
use as excitation dynamics a SIRS dynamics. A suscep- 
tible agent gets the disease by keeping physical contact 
with an infected one for certain time. The transmission 
of the infection (=excitation) is modeled by a Poisso- 
nian process whose characteristic time is tt- Once in- 
fected, an agent carries the infection for a certain period 
of time, after which gets recovered and remains immune 
to the disease for another period of time, to finally be- 
come susceptible of being infected. The transitions from 
infected to recovered, and from recovered to susceptible 
are modeled by Poisson processes characterized by the 
characteristic times t/ and r^j, respectively. 

Low density. — At low density the agents do not form 
large clusters and we can reduce the dynamics to a much 
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FIG. 2: Populations of agents at the steady state vs. linear 
system size L. Simulations performed with the parameters 
given in fl3] and v = 0.1. Full black lines corresponds to Eqs. 
dS]), and ^ while the open symbols to S"'^ , ift^ and 
R?t^ ■ The vertical line corresponds to Lc. 



simpler one in which we consider only binary collisions. 
In this gas-like phase we describe the time evolution of 
the populations of susceptible, infected, and recovered 
agents by the following mean-field equations: 
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(3) 
(4) 
(5) 

where S, I and R are defined as S = Ns/N, I ^ Nj/N 
and R = Nji/N, with N being the total number of agents 
in the system, and Ns, Nj, and Nr being the number of 
susceptible, infected and recovered agents, respectively. 
Here, ctq represents the scattering cross section of the 
agents, defined as do — 4r, p is the density (p = N/V), 
and uj denotes the mean duration of the collision event. 
The exact dependency of lo with v depends on the par- 
ticular agent-agent interaction rule. Since in general, for 
V — > 00, UJ has to vanish, while for v 0, tu 00, and 
considering to the Buckingham tt theorem plj . we con- 
clude that UJ takes the functional form uj (v) = K (r/v)^ , 
where ^ is a positive constant and K a dimensional func- 
tion of the interaction rule parameters. We choose values 
of K and ^ that best fit simulation data of uj vs. v [l^. 
In consequence, (l — e"'^^'^^) represents the probability 
per contact event (of mean duration uj) of transmitting 
the disease, while vaop in Eqs. ^ and ^ represents the 
agent-agent contact rate (CR). 

From Eqs. ©-([n]), we obtain the two steady states of 
the system: i) 5 = 1, / = 0, and i? = 0, corresponding 
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to the epidemic extinction, and ii) 

Sst = Rq ^ 



1st = (l-Ss 

Rst = (1 ^ Sst) 



Tl 



(6) 
(7) 

(8) 



corresponding to the endemic state, where Rq = 
(l — g^-^/'^T^ ^^^p^^ denotes the basic reproductive num- 
ber. Recall that Rq > I corresponds to the endemic state, 
while i?o < 1 to epidemic extinction. It is instructive to 
represent Rq as Rq — X{k), where (k) denotes the mean 
number of individuals in contact with an infected agent 
per time unit, i.e., (k) = vctqpti, while A refers to the 
probability for a susceptible individual to contract the 
disease during a contact event with an infected one, i.e., 
A = (l — e^'^/'^^). Notice that though A and (k) are 
typically independent model parameters, for the current 
mobile agent-based model, they are coupled. This cou- 
pling has proved to be essential to understand the disease 
spreading dynamics. 

Fig. [2] shows a comparison of Eqs. ©-([H]) and agent- 
based simulations. For large values of L, there is a good 
agreement between simulations and theory. For small 
values of L, spatial correlations become important and 
multiple collisions and shielding effects take place. How- 
ever, since the spatial dynamics is faster than the exci- 
tation dynamics, the mean-field approach still approxi- 
mates the simulation data as observed in Fig. [21 Clearly, 
the smaller L, the worse the approximation becomes. No- 
tice that the critical system size Lc that separates the 
endemic from the epidemic extinction regime takes the 

1/2 

form: Lc — \N (l — e^'^/'^'^) vaQTi\ (dashed vertical 
line in Fig. [2]) . 

High density. — At high enough densities, agents form 
percolating clusters and agent mobility is strongly re- 
duced, though exchange of neighbors still takes place. 
The maximum conceivable density, corresponds to 
the one associated to the maximum packing fraction 
of disks, TjM = Pmcl, where a is the area of a sin- 
gle agent. rjM is 0.91 and corresponds to an hexago- 
nal lattice. At densities close to this limit, the system 
can be reduced to a regular lattice of agents through 
which a SIRS dynamics evolves. Since we assume that 
locally there is still some mixing of particles, we ne- 
glect correlations between nodes, and express the tem- 
poral evolution of agent-populations in terms of a sim- 
ple mean-field which yields the following endemic states: 
S^,^ - tt/{tik), = (l-S,t){Ti/iTj+Tn)), and 
R^t^ = (1 — Sst) {tr/ (t/ + tr)), where k denotes the 
mean number of neighbors. These values are shown in 
Fig. [2] as open symbols and represent a rough estimate 
for the system behavior at densities close to pM- Bet- 
ter estimates can be obtained by considering spacial cor- 
relations through the two-point probability function as 



o 



^ = 1.3 
f = 1 
^=0.7 



R 

MF Theory 



FIG. 3: Populations of agents at the steady state vs. agent 
active speed v. Simulations performed with the parameters 
given in [l^ and L = 150. Full black lines corresponds to 
Eqs. ©-([H])). The vertical lines correspond to Vm and Vc, 
while dashed, dotted, and dash-dotted curves correspond to 
different values of 5, see text. 



indicated in [15|. Notice that k is constrained to the 



structure of the lattice and that S^^ > 0. 

The role of the active speed. — The active speed of the 
agents u is a control parameter which directly regulates 
the degree of interaction among agents. Particularly, CR 
is directly proportional to v. The values of the epidemic 
steady state, as indicated by Eqs. (I6|)-(l8]) and confirmed 
by agent-based simulations, exhibit a non-trivial behav- 
ior with V (see Fig. [S]). For low enough density, we expect 
to observe the epidemic extinction state in the limit of v 
going to zero. This implies that to observe the endemic 
state, V has to be larger than a threshold value Vm- From 
Eq. ([6]) we can estimate the minimum speed required by 
the disease to survive. This is set by the condition Sst = 1 
from which we get the following transcendental equation 
for the minimum speed v^'- 



1 



1 



(9) 



By solving numerically Eq. ^ for the parameters used 
in Fig. [21 we obtain Vm — 0.027 (dashed vertical line). 
Note that the expression for Rq predicts a crossover at 
Lo{vc) = tt (dash-dotted vertical line). 

As shown in Fig. [31 there are three regimes with the 
agent active speed v. For v < Vm, the mean time be- 
tween subsequent collisions is larger than the infection 
period r/, and so the disease can not spread and gets ex- 
tinguished. On the contrary, for < v < Vc, the mean 
time between subsequent collisions is smaller than r/ and 
the system ends up in an endemic state. In this range 
of V, the mean collision time is such that uj{v) ^ tt, 
and thus the probability of transmitting the disease in a 
collision event is very high. As a consequence, the epi- 
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demic size increases with the CR leading to Sst ^ v~^. 
Notice that in the Hmiting case of tt 0, this regime 
is vahd for all v > u,„. However, any realistic disease 
has some non- vanishing characteristic transmission time, 
and consequently, a third regime, ignored so far, emerges 
for V > Vc, i.e., for larger values of CR. In this case, 
to{v) < Tt and the mean collision time is not enough 
to assure the disease transmission in each collision event. 
The interplay between CR and the transmission probabil- 
ity per collision event determines the third regime which, 
contrary to the second one, is model dependent and such 
that Sst ^ v^^^- In consequence, there are three possible 
behaviors depending on the interaction rule, resp. ^. If 
^ = 1, in the limit of u — > cxd the epidemic size reaches an 
asymptotic value Soo ~ tt/ [KraoTjp] > 0. For ^ > 1, 
the disease spreading exhibits a non-monotonic behavior 
with V and in the limit o^v—^oo the systems reaches the 
epidemic extinction regime. Thus, an increase in w, resp. 
CR, leads to a reduction of the epidemic size. Finally, for 
< ^ < 1 and in the limit of u — > oo, Sst goes to zero, 
which indicates that the epidemic size at low density and 
for large enough values of v can become larger that the 
bounded epidemic size corresponding to S^^ . Note that 
simulations performed with the specific agent-agent in- 
teraction rule given by Eq. ([2]) fall into this case, with 
C = 0.91. 

Concluding remarks. — We have shown that in mo- 
bile agent systems the coupling between mean exposi- 
tion time and CR is crucial to understand the excita- 
tion dynamics, which generically exhibits three regimes 
with CR: excitation extinction for low CR, an excited 
regime where the excitation is such that the number of 
quiescent agents is inversely proportional to CR, and a 
third regime, for high CR, where the number of quiescent 
agents scales with an exponent ^ — 1, with ^ being the 
scaling exponent of the mean exposition time with CR. 
The novel third regime is clearly model dependent and 
opens, in the context of epidemic spreading, the counter- 
intuitive possibility of " curing" by increasing CR, which 
can be achieved in simulations by modifying the interac- 
tion rule among agents such that ^ > 1. 

At the experimental level, the introduced mobile agent 
model might help to understand spatial distribution (and 
levels) of gene expression induced by cell-cell signaling as 
observed in some bacterial colonies 
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It also might 

shed some light on disease spreading on bacteria l23l . [2 3 
and also on competing strain bacteria experiments [25| . 
Clearly, the current agent model is not suitable to de- 
scribe disease spreading on structured populations, i.e., 
when the disease moves across metapopulations, as in 
humans. Further improvements to account for struc- 
tured communities might include preferred interacting 
partners, which can be modeled by a long-range poten- 
tial acting among members of a community, as well as by 
dividing the space into several connected compartments. 
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